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ABSTRACT 

We present an analytic derivation of the temporal dependence of the perpen- 
dicular transport coefficient of charged particles in magnetostatic turbulence, for 
times smaller than the time needed to charged particles to travel the turbulence 
correlation length. This time window is left unexplored in most transport mod- 
els. In our analysis all magnetic scales are taken to be much larger than the 
particle gyroradius, so that perpendicular transport is assumed to be dominated 
by the guiding center motion. Particle drift from the local magnetic field lines 
and magnetic field lines random walk arc evaluated separately for slab and 3D 
isotropic turbulence. Contributions of wavelength scales shorter and longer than 
the turbulence coherence length are compared. In contrast to slab case, particles 
in 3D isotropic turbulence unexpectedly diffuse from local magnetic field lines; 
this result questions the common assumption that particle magnetization is inde- 
pendent on turbulence geometry. Extensions of this model will allow for a study 
of solar wind anisotropies. 

Subject headings: 
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Introduction 



The behaviour of individual fast charged particles in magnetic turbulence is relevant to 



a number of problems in plasma astrophysics, from the solar wind (e.g. Bruno & Carbone 



2005) to interstellar medium (e.g. Elmegreen & Scalo 2004) and cosmic rays at highest 



energy (e.g. Fraschetti 2008). However, in contrast to cosmic rays with energies beyond the 
GeV scale, a thorough understanding of the particle transport properties can be attained 
only in interplanetary space, where in situ measurements of both magnetic turbulence 



energy spectrum and particles energy are possible. Diffusion theory (e.g. Jokipii 1966), as a 



main tool to study charged particle propagation in magnetic turbulence, yields a statistical 



description of a population of particles and relies on the approximation (Jokipii 1972) that 



a characteristic time T exists much larger than the correlation time tc of the magnetic 
field fluctuations (as seen by the particle) but also much smaller than the time-scale of 
both the variation of these fluctuations and of the average distribution function. The 
Vlasov-Boltzmann equation for the charged particles phase-space distribution function can 
be therefore considerably simplified to terms of the second order moments of the magnetic 
field fluctuations. In this scenario higher-order moments are not necessary to determine 
the particles motion as the process is markovian; diffusion is governed by the central limit 



theorem ( Chandrasekhar 1943). 



Observational constraints posed by heliospheric environment on perpendicular diffusion 



across the average magnetic field involving, e.g. jovian electrons (Chenette et al. 1977) 



have not yet been included in a first-principles unified theoretical picture. Perpendicular 
diffusion occurring in the ecliptic plane is invoked as a plausible explanation of the time 



delays in solar energetic particle events detected by Helios (Wibberenz & Cane 2006). A 



more remarkable longitudinal separation in the combined electron observations by Stereo 
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A/B and SOHO from the January 17, 2010 event [] suggests a strong diffusion perpendicular 
to the mean magnetic field. The access of high energy particles to high-latitude heliospheric 



regions observed by Ulysses (Malandraki et al. 2009) is dominated by particle propagation 



along the mean magnetic field lines, although cross-field diffusion cannot be excluded. 

Perpendicular transport in a magnetic field depending on two or fewer space coordinates 



originates only from the meandering of the magnetic field lines (Jokipii et al. 1993) whereas 



in an arbitrary three-dimensional turbulence also emerges as a general property of the 



particle motion (Giacalone & Jokipii 1994). Kinetic approach has been applied to 



perpendicular scattering of strongly magnetized charged particles by using a model for 



the collision integral (Chuvilgin & Ptuskin 1993). Recent numerical simulations (Minnie 



et al. 2009) investigated the common assumption that the charged-particle gyrocenter 



follows the magnetic field lines: the approaching of the guiding center cross-field motion 
to the transverse field line random walk for various parallel mean free paths is studied 
for a solar-wind like turbulence. However we notice that the assumption that guiding 
center follows the magnetic field lines does not result directly from the equation of motion; 
therefore it may be realized only for particular turbulence models. 

An approximate diffusive perpendicular transport model based on the guiding center 



motion (Non-Linear Guiding Center, NLGC) has been put forward in Matthaeus et al 



(2003), which provided a method to compute magnetic fluctuations along perturbed particle 



trajectory. However, the NLGC's assumption that the probability density of perpendicular 
displacement is diffusive at all times is a limitation of this model. Subdiffusive nature of 
perpendicular transport in a slab turbulence was not recovered, contrary to the expectation 



from the conservation of canonical momentum of ignorable coordinate (Jokipii et al. 



1993) and to the findings from test particle numerical simulations in turbulence having a 



^Available at URL http : / /w'W'w2.physik.uni — kiel.de/ stereo/ downloads / sept-electron.events.pdf 
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dominant slab component (Qin et al. 2002). NLGC has been extended to early phase of 



perpendicular scattering (le Roux et al. 2010), where the probability distribution function 



of the waiting-time between two scatterings decays slower than the exponential, therefore 
including memory effects (non markovian process). For a review of other improvements of 



NLGC see references in le Roux et al. (2010). 



In this paper, we explore the transition to the diffusion regime in a magnetostatic 



turbulence by using the first-order orbit theory as proposed by Rossi & Olbert (1970), which 



is based on two assumptions: 1) the particle gyroradius is much smaller than any variation 
length scale of magnetic field and 2) the turbulent magnetic energy is much smaller than 
the average magnetic field energy. This allows us to disentangle the cross-field particle 
motion into two separate components: field lines meandering and gradient/curvature drift 
from the local field line. In the present paper, the drift is meant to be the transverse 



gyroperiod- averaged motion of the guiding center away from the local field line (Rossi & 



Olbert|[l970 ). Field- line meandering has been recognized in early solar turbulence studies to 



be the main contribution to perpendicular motion to average field direction (Jokipii 1966 



Jokipii & Parker 1968); the diffusive nature of the field line spread has been related to the 



turbulence power spectrum power law index at low wavenumber (Ragot 1999). On the 



other hand, individual particle gradient /curvature drift on scales smaller than correlation 



scale does not seem to have been object of theoretical investigation (see however Schlickeiser 



& Jenko (2010)), on the basis of the common belief that particle magnetization does not 



depend on the particular turbulence geometry, in contrast to recent numerical findings 



(Tautz & Scalchi 2010) Therefore, the knowledge of particle trajectory along and across 



local field lines has remained completely undetermined. In this paper we shed light to this 
distinction analyzing cases propedeutically relevant to the solar wind propagating cosmic 
rays and astrophysical blast waves of supernova remnants. 
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2. Transverse guiding center drift 

We consider a process of propagation of a charged particle in magnetic turbulence 
which is statistically homogeneous in time, i.e., the velocity correlation depends only on 
time difference along the orbit. The instantaneous mean square displacement along the 
space coordinate x after a time At = t — t' for a particle propagating in an arbitrary 



Taylor 


1922 


Green 


1951 


Kubo 


1957) 



Jo 



where t' is an arbitrary initial time, ^ the time lag and the ensemble average (..) is meant 
to be the average over a population of particles and over an ensemble of turbulence 
realizations. We notice that Eq.Q applies for any value of At. We define 



dxAt) ^ --{(AxY) 



(2) 



The standard perpendicular coefficient of diffusion can then be deffned as 



Kfr 



lim dxxit). 



(3) 



We here investigate the transverse motion of a low-rigidity particle for a time smaller than 
tc in general 3D magnetostatic turbulence. The diffusion approximation therefore may not 
be valid. 

We consider a spatially homogeneous, ffuctuating, time-independent magnetic field. 
The amplitude of the fiuctuation (SB) is assumed to be much smaller than the average 
field magnitude (Bq). We represent such a magnetic field as B(x) = Bp + 5B(x), with an 
average component Bq = BqBz and (5B(x)) = and 6B{x.)/Bq ^ 1. This approximation 
is known to be valid in several turbulent media, as the solar wind, where the propagation 
of the magnetic fiuctuation is much smaller than the velocity of the bulk ionized fiuid. We 



will make use of the first-order orbit theory (Rossi & Olbert 1970): the particle gyroradius 
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r„ is much smaller than the length-scale of any magnetic field variation: 



To <^ mm 



(4) 



where Bi is the z-th component of the perturbed field B(x). No further assumption is 
made on the spatial dependence of 5B or geometry. In this approximation, we consider 
the guiding center motion. In a spatially varying magnetic field, the guiding center may 
significantly drift from the average field direction due to the action of the field gradient 
on the particle magnetic moment. We therefore consider non-zero gradient and curvature 
drifts. We estimate that drift and resulting displacement after a time shorter than the 
correlation time. The guiding position X(t) = (X, Y, Z) for a particle of mass m, charge Ze 
and momentum p having coordinate x(t) = {x, y, z) is described, in c.g.s. units, by 

c B X p(t) 



X(t) = x(t) 



(5) 



Ze 52 

If the scales of magnetic fiuctuation are much larger than the gyroradius r^, the guiding 
center motion defined in Eq. (|5]) has the role of effective gyroperiod-averaged motion. 
Therefore, Eq. (|5]) can well describe the motion perpendicular to the local magnetic field. 
In the case of "finite Larmor radius", the gyroradius only represents the typical scale of 
particle motion and Eq. ([s]) provides the instantaneous guiding center position whereas 
gyroperiod-average becomes meaningless. In the magnetostatic field described above, the 
guiding center velocity transverse to the field B(x) is given at the first order in 6B(x.)/Bq 



by the gyroperiod average (|Rossi & Olbert||1970|) 
Vlit) 



vpc 
ZeB^ 



1 + 



B X VB + ix'B{V X B)_ 



(6) 



where a is the particle pitch angle and /i = cosa. Eq. ^ gives the first order most general 
expression of the guiding center velocity orthogonal to the local magnetic field direction 



(Balescu 1988). Here the variation of a is assumed to be negligible over a gyroperiod. 



Being \j_{t) a gyroperiod average, magnetic field can be computed at the guiding center 
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position during that gyroperiod. In contrast to Matthaeus et al. (2003), the transverse 
motion of the guiding center from the field hne is not parametrized in the present paper 
through some constants to be inferred from numerical simulations, but described directly 
from the equation of motion of the guiding center. The finite-time average square transverse 
displacement of the particle from the direction of local B due to drift doit) can then be 
written in this approximation using the Eq.(|2|: 



Jo 



(7) 



where i stands for any transverse coordinate, X or Y. The average square displacement is 
computed from the following expression, to the lowest order in 6B/Bq, 



f vpc \ ' 



da 



djSB^ + fi^dsSBj 



[x(t')] X 



-dj6B3 + fi^dsSB, 



where = (1, 2) or (2, 1) and the fields are evaluated at the perturbed particle position 
x(t). The average square of the displacement du-^it) does not depend on the sign of 
the electric charge Ze, at variance from the drift velocity in Eq. We notice that 
in the first-order orbit approximation the particular case of 2D turbulence defined by 
SB{x,y) = {6Bx,6By,0) provides a zero transverse velocity drift; this is because, on the 
right hand side of Eq. (|8|, this form of turbulence has SB^ = and also 5Bj does not 
depend on the z coordinate. Thus, this analytic method cannot be applied to the composite 



slab/2D solar wind model of Bieber et al. (1996), a very useful but empirical description 
of the MHD-scale turbulence in slow solar wind. Moreover, the slab/2D model could 
be incomplete as the non-wave ("2D") turbulence might have an additional component 
along Bq. This implies that, due to the sub-diffusive nature of the perpendicular particle 
transport in slab turbulence, drifts from local field-line found in numerical simulations for 
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composite model (Minnie et al. 2009) are second-order contributions. Different anisotropies 
may be compatible with large-scales solar wind observations; in this paper we indicate a 
possible alternative method. 

We may simplify the derivation by using the Fourier representation of (5B(x): 

/oo 
£k5B{k)e'^<''^ (9) 
-oo 

where 3?(-) stands for the real part and x(t) is the particle position at time t. Therefore the 
average displacement in Eq.(|8]) contains terms of type 

/oo 
d'^k5Bj{k){iki)e'^-''^'^ (10) 
-oo 



with l,j = 1,2,3. We compute the particle position in Eq.(lO) along the local 
magnetic field: x(t) = xo(t) + :s.MFL{z{t)), where the unperturbed particle orbit is 
xo(t) = (fsin0A/l — yu^/f^, —vcos(j)^yi — yU^/fi, t>||t); here v\\ is the unperturbed particle 
velocity along z, (j) is the particle azimuth angle in the plane orthogonal to Bq and 
Vt = ZcBq / [ni'-yc) the particle gyrofrequency in the background field containing the Lorentz 
factor 7; ^MFL{,z{t)) = {xmfljVmfl, z{t)) is the offset in the plane orthogonal to Bq due 
to the magnetic field line random walk (MFLRW) at z = z{t). The assumption of ballistic 
motion along Bq, i.e., z = fyt, relies on the choice 6B ^ Bq] at times smaller than the 
correlation time of the perpendicular fluctuation, a fortiori we cannot assume parallel 
diffusion. At the small length-scales considered here, parallel and perpendicular motions can 
be disentangled and any non-markovian parallel motion, e.g., memory effect of a particle 
tracing back its trajectory, is not expected to interfere with the perpendicular transport, in 
contrast to the case of compound diffusion. We can write 

The magnetic field lines (MFL) are defined by dyiMFL x B = 0. This implies that 
a finite distance A^s^^jfl in the ballistic approximation is a first order term in 5B: 



Axmfl — ('^B/i?o)f lit. Therefore, the exponential e*'^'''*^^^*^^*^*^-* contribute only at zero 



order in Eq. (11) and the fuctuation in Eq. (10) can be computed along the unperturbed 
trajectory: e''^'''^*^ ~ g«k-xo{t)^ which is equivalent to the quasi-linear approximation. To first 



order in Eq.(lO) we can replace the exponential as 

gik-x(t) ^ gi(VKsin(i/;-(/))+fc||t)||t) 



(12) 



where = tg ^{ky/kx), k\i = kz, k± = ^Jk1 + ky and W = k±v^yi — jj'^/Q = k±rg. 

By using the Bessel function identities (see [Abramowitz &: Stegun (1964), Eq. (9.1.41)) 

oo 



Eq.(ll) is rewritten as (see also Schlickeiser (2002), Sect. 12.2.1): 

oo 

gik-xo(t) _ j^^|y^giA:||i)||t+in{?/i-</)+™) 



(14) 



The magnetic fluctuation space derivative in Eq.(lO) can be then written in the following 

way 

di6Bj{^) = ^ J2 d^k5Bj{k){iki)J^{W) X 

„ ^ — oo 



with /, j = 1, 2, 3. The typical term in Eq.([8| is of type 



(15) 



here F^jj?) represents various /i factors resulting from the expansion of Eq.dsl). Using 



(16) 



Eq.(15), we obtain for Eq.(|16|) 

3? 



F(/i2)x 



oo oo 



d^iY^ Yl / / d^k'5Br{k) X 

n=-oo m=-oo "^-'^ -^-"^ 

Uw)itk)5B;{k')j:xw)i-tk'^) X 



(17) 
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We assume the standard inertial range magnetic turbulence power spectrum which is 
uncorrelated at different wavenumber vectors: 

{SBr{k)SB*(k')) = 6{k - k')Pr,{k) . (18) 



Thus Eq.(17) reduces to 

nf^') E J d'kR{k,t)Pr,{k)kk,4{w) 

whose time-dependence is entirely contained in 

ft „-i(A;||D||+nf2)t i 

i?(k,t) = / rfee-^('=ii^ii+"^)« = ^ . 

Jo -t{kiivi\ + nil) 

Since '^Jn{W) = 0, where stands for imaginary part, we may consider ^R(k,t): 

II + nil 

The orthogonal scale l/k± can be estimated as \Bi/djBi\, thus Eq.Q states 

W ~ k±v/n < 1 . 



(19) 



(20) 



(21) 



(22) 



For W <^1, it holds Jo{W) > Jn{W) for n > 1; moreover, ^R{k,t) ~ 1/n for large n. We 



may therefore approximate the sum in Eq.(19) as its term with n = 0. Therefore Eq.(19) 



yields, using Eq.(21), four terms of type 

2 



/ vpc 



oo 

F(/i2) / d^kPr,{k)kik, 



sin[A;||f ||t] 

ki\v\\ 



(23) 



with indexes {r,q,l,p) = (3,3,2,2), (3,2,2,3), (2,3,3,2), (2,2,3,3) for doxx 



{r,q,l,p) = (3,3,1,1), (3,1,1,3), (1,3,3,1), (1,1,3,3) for dDyy- -^l- (23) represents the 



general term contributing to the first-order transverse drift coefficient of a particle in a 
static first-order perturbed magnetic field. 
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3. Magnetic-field-line random walk 

In the present section we compute the contribution to the time-dependent particle 
transverse transport due to MFLRW. If the correlation function of the magnetic fluctuation 
is homogeneous in space, the mean square displacement of the MFL orthogonal to the z-axis 
can be defined, in analogy to Eq.([2]), as 



^ rdz'{6B4^iz')]6B,[^m)- (24) 



We compute the magnetic turbulence 5B(x) in Eq.(24) along the unperturbed trajectory of 
a particle travelling with zero pitch-angle, as in Eqs.([9, 12). The motion along the average 
field is then ballistic, i.e. z = v\\t. In these approximations the transverse displacement of a 
MFL corresponding to a distance v\\t along Bq travelled by a small rigidity particle can be 
written as 



°° sml 



dMFLit) = I d'kPr.ikf-^^^^^ . (25) 



Equation (25) is in agreement with Eq. (17) of Shalchi (2005) derived for pure slab 
turbulence. The MFL coefficient diffusion describing the random walk of the field lines can 
then be defined as 

hmfl = lim dMFiit). (26) 

t— s>oo 

In our approach, MFL diffusion cannot be assumed because travelled distance z smaller 
than parallel correlation lengths is considered; nevertheless, a simple ballistic motion along 
the z-axis allows to recover the standard result of MFL perpendicular diffusion in QLT. The 
discussion of the previous two section implies that the instantaneous coefficient diffusion 
perpendicular to the average magnetic field Bq is given, in the presence of weak turbulence 
and neglecting parallel scattering, by two contributions: the random walk of the field line 



- 13 - 



and the guiding center drift from the field hne: 



d{t) = doit) + v\idM Flit), 

K = kd + v\\kmfl = ^i'^ldoit) +v\\dMFL{t)]- (27) 



t—>-oo 



In the next section the previous results are applied to the slab and 3D isotropic turbulences. 

4. The turbulence power spectrum 

In this section we apply the approach developed in previous sections to derive the 
instantaneous transverse particle transport coefficients both of guiding center drifting from 



local MFL (see Eq.(23)) and of the MFL from the average field direction in Eq.(25) by using 



the coherence length of the turbulence to disentangle small from large scale contributions to 



perpendicular diffusion. We will consider two cases: 1) slab turbulence, introduced (Jokipii 



1966) to represent the static limit of the solar wind magnetic fluctuations and extensively 
studied with Monte Carlo numerical simulations; we will compare the result with the QLT 
limit; 2) 3D isotropic turbulence, idealized case likely to provide an unperturbed model for 
anisotropics observed in the solar wind. 

4.1. Slab 

We consider the slab turbulence, static limit of transverse and longitudinal- 
propagating Alfven waves: (5B = 5B(z) and 5B(x) ■ e^ = 0. In this case, from Eq.(|8]), 
^Dxx^^^ ~ ^Dyy^^^ ~ ^hi^)- The turbulence wave number is aligned to the average magnetic 
fleld, thus we adopt the following form of the power spectrum: Prg(k) = G{k\\){6{k±)/k±)6rq 
with r,q = 1,2, and P3i(k) = with i = 1,2,3. The ID spectrum is assumed to be of 
Kolmogorov type. Observations of electron-density fluctuations inferred from scintillation 
measurements exhibit a Kolmogorov power law, with index approximately equal to 5/3, 



over 5 orders of magnitude (Armstrong et al. 1995). Several other observations of magnetic 



turbulent media, from earth's magnetosphere to galaxy clusters, validate the Kolmogorov 
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power spectrum up to a range of 12 orders of magnitude. We mention that solar wind 
observations show that at scales smaller than the ion thermal gyroradius (~ 10^ cm around 
the earth), much smaller than the scales considered in this paper, the magnetic turbulence 



spectrum deviates from the Kolmogorov, having an index of —2.12 (Bale et al. 2005). 



At length-scales larger than the coherence length the measured interplanetary magnetic 



turbulence is well described by a flattening power spectrum (Hedgecock 1975 Bieber et al 



1994). On the other hand, a consistent comparison with the quasi-linear limit requires the 
power spectrum to be defined at scales larger than coherence length, i.e. for A;|| < A;™", up 
to the physical scale of the system 2Ti/k^^^ \ we will adopt here a simplified form: 

G'(fcii) = (28) 
[ G'jj(A;™")-'? if /tf < k\\ < fc™" , 

where /c™"^ corresponds to the scale where the dissipation rate of the turbulence overcomes 
the energy cascade rate. The choice of a constant power spectrum at large scales instead 
of a function smoothly connected to the inertial range already used in the literature is 

^0 



merely dictated by easier mathematical tractability. Here g = 5/3 and the constant G?, is 



determined from the normalization 

/oo 
£k{Pn + P22 + P33) (29) 
-00 

implying, using cylindrical coordinate {d^k = dk\\k^dk^dip), 



(30) 

II A t l„min\\—q ' ^ ' 



^:'nq{kl 



with the assumption k\ ^ A;?^™ ^ A;, 



max 
II ■ 



We consider first the transverse drift in Eq.([23|). We average -F(/i^) over an isotropic 
pitch angle distribution. Using cylindrical coordinate [d?k = dk\i^k±dk±dip) we have 

-y^^ dkfyGik,)^— (31) 
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In units of the Bohm coefficient diffusion [kb = (l/3)rgf) and approximating r^/fn ^ ^, 
we obtain 



dUt) _ 3 (SBVq-l 



(32) 



where we defined 



-^(2 - sint/ii - yn cosyy 



(33) 



where the time-dependence is contained in the new variable = k^'^"'v\\t ~ k]^'^'^rgVLt (and 
= ||t — k^\^rgVLt) and we used 



/(a, m) 



^sinydy 



i/2[e-'^''T{a, lu) - e*5"r(a, -m)] 



(34) 



where r{a,z) is the incomplete gamma function (see [Gradshteyn fc Ryzhik| ( [l973| , 
Eq.(3.761.2)). The time evolution of the drift coefficient d'}^ is depicted in Figjl] The 



ffist term in Eq. (33), corresponding to scales smaller than coherence scale 27i/k', 



{k > fcji"*"), dominates over the second term, corresponding to scales larger than 27[/k'!^" 
{k < A;™"). The diffusive behaviour can be found by using the approximation of r{a,z) 
for \z\ = y^ = /c™"f||t ^ 1 (see 



Gradshteyn & Ryzhik (1973), Eq.(8.354.2)), because 



^ 1; since we assume a weak magnetic fluctuation, it is reasonable to assume that the 
perpendicular diffusion time-scale is shorter than the parallel scattering time-scale, i.e.. 



1/A;|l^*"t'|j , or in other terms the diffusion limit is the dominant term in Eq. (33) for large t 



and t < 1/A;™"t;||: T{a,z) ~ T{a) - z^/a; therefore J(2 - g,y||) ~ sin(g7r/2)r(2 - q) (dashed 
line in FigjT|. In diffusive regime, the second term in Eq.(33), representing the large scales 
(/ > 271/^™"" or k < fc™"), does not significantly contribute to the particle drift, as it is 
manifest in Fig. [T] We find that for slab turbulence, transverse particle drift coefficient 
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from local MFL is given by 

dm 3 f6By q-lsm{qn/2)T{2-q) 
KB 20\BoJ q (A;™%)i-9(fit)2-9 ' ^ > 

thus subdiffusive with behaviour ~ t^^'^^'^^ (depicted as the dashed line in Fig. [l 

Transverse subdiffusion has also been found by considering time-scales longer than the 



parallel scattering time and therefore allowing parallel scattering in Kota & Jokipii (2000). 



However, in that case particles are assumed to propagate back and forth along the MFL 



and to be tied to the MFL. We notice that the drift-coefficient time evolution in Eq. (35) 



confirms that charged-particles in a turbulence depending on less than 3 space coordinates 



remain confined within a gyroradius from the local field line (Jokipii et al. 1993; Jones et 



al 



1998). The time-integration of Eq. (35) up to t = Ly/i; ~ 2tc /{kl^^'^rgQ), gives in case of 
weak turbulence {6B ^ Bq) the condition (Ax^) <^ r|. We notice that the time-integral of 
k|) ~ t^^^, which provides ((Ax)^) ~ t^^^, is an increasing function of time for any observed 
physical value of q; however, as shown above, this result does not contradict the theorem 
of reduced dimensionality. In summary, the present result has been obtained under three 
assumptions: 1) ballistic motion in the z coordinate {z = v\\t)] 2) average displacement 
transverse to the local field B due to first-order drift; 3) Kolmogorov power spectrum for 



magnetic fiuctuations. Equations (32, 33) represent the average transverse displacement 
computed in the first-order orbit approximation at any time smaller than the parallel 
scattering time-scale, so that the approximation of ballistic motion parallel to the mean 
magnetic field holds. 



From Eq. (25), the MFLRW in units of magnetic coherence length Ly = 27r/A;™"' is 
given by 

d-MFLiW = (^)' Sr^^^"'^^ ^^^^ 

where we defined 

i7(y|,, g) = {kf^r,ntyi{-q, + Si{yr^) , (37) 
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where Si(a;) is the Sine integral function. The first term in Eq. (37), corresponding to 
scales smaller than coherence scale 27r/fc?j"*" {k > kJP'^^), is dominated by the second term, 



corresponding to scales larger than 271 /kV^^^ {k < A;|^™). From Eq.s (36, 37), the MFLRW 



diffusion coefficient is given by = {^B / BQY'n{q — 1) / {Aqk'^'^'^) . In Fig. [2| the i^^Mp^ is 



shown to recover the quasi-linear limit and is dominated by large wavelengths, given by the 



Si(a;) term in Eq. (37): Dmfl = tt G(k\\ = 0)/Bo = i^%,ifl^ where the quasi-linear limit is 



expressed, as known, as power spectrum at zero parallel wavenumber. 



Equations (36, 37) provide the MFLRW for distances Az along Bq shorter than Ln 



The guiding center perpendicular scattering in a slab turbulence is described as a series of 
bumps in the transverse drift which are asymptotically suppressed confining the transverse 
motion to follow the MFL meandering, as also found in low-rigidity test particle numerical 



simulations (Qin et al. 2002). Therefore, we confirm that slab transverse transport is due 



to the meandering of MFLs but we also model the transport across the MFL for first-order 



magnetic fluctuations not considered in previous treatments (Kota & Jokipii 2000). For 



the slab turbulence, the transverse particle diffusion can be disentangled in two energetic 
contributions: drift coefficient is dominated by length scales smaller than coherence length 
whereas the MFLRW is dominated by length scales larger than coherence length. 



4.2. Isotropic 



We consider 3D isotropic turbulence, in which the turbulence 5B depends on all three 

trui 

kik„ 



space coordinates. We adopt the following power spectrum (Batchelor 1970): 

G{k) 



with 



G{k) 



G^k-i 
Gt\k^- 



if kr, 



k^ 



<k <kr, 



if kn < k < kr^ 



(38) 



(39) 
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where corresponds to the scale where the dissipation rates of the turbulence overcomes 
the energy cascade rate, the coherence length is given by L = 27r/A;™", and the physical 
scale of the system by 2Ti/k^ (about the spectrum at large scales, see the discussion in 



Sect, 4.1). The constant Go is fixed by normalization: 



C„.&zi). (40) 

assuming ^ ^f^m ^ ^max_ spherical coordinate for the wavenumber 

k = A;(sin 6* cos sin6' sin cos 6') and 3D-turbulence: 5B(x) = {6Bx,SBy,6Bz). We 



compute first (i}) (t). In reference to Eq. (23 ), the non- zero terms are {r,q,l,p) = (3,3,2,2), 



{r,q,l,p) = (2,3,2,3) and {r,q,l,p) = (2,2,3,3). Combining the non-zero terms and 
averaging over a pitch-angle isotropic distribution, this gives 



Kb 4:\BoJ q 



d\){t) 3 f6By q~l 



1 ikmin^ 



a) 



X 



+ 



yo 



ym Um 



(41) 



with Um = kminV\\t ~ k„iinrgVLt (and Uq ~ k^rgVLt)] here, the term integrated over scales 
smaller than coherence scale 27r /kmin (or k > kmin) is recast as 

Fi{y,q) = —- [cosy + ySi{y)) + 

15 / — g 5 g 

Uy-^-'^smy 3q^ - Uq + 18 - 

- y^T^ + ' 15g(2-g) '^'-^'^^ 

14 1 + (7 , , , 

and the term integrated over scales larger than coherence scale 27i/kmin (or k < kmin) is 
recast as 

15 5 

21 — hy"^ 7 cosy 
-|- smy (43) 
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where we approximated Tg/v\\ c^i^Vt ^ as in Eq. (32) and we used 

I{a,u) = I y°'~^cosydy 



l/2[e-*5'^r(a, iu) + e'5"r(a, -iu)] , 



(44) 



from Gradshteyn & Ryzhik (1973), Eq.(3.761.7). The instantaneous transverse coefficient 



diffusion in Eq.(41) is represented in Figj3j The dominant term in the diffusive hmit, with 
the condition ym = kminV\\t < 1, is given by 



K 



D 



(45) 



Kb AO \ Bq J q — 2 " Kb 

and represented in Figj3] In contrast to the slab, the particle drift from the MFL does 
not depend only on the power spectrum at length-scales smaller than Ly. As for the slab, 
transverse diffusion is axisymmetric: (^^^^^(t) = ^^^^^(t) = ^^^(t). Statistically, charged 
particle motion is not tied to local MFL in 3D isotropic turbulence. Theorem on reduced 



dimensionality turbulence in Jokipii et al. (1993) and Jones et al. (1998) allows charged 



particle to be magnetized to local MFL within gyroradius scale only in turbulence depending 
on a reduced number of space coordinates. Any three-dimensional extension could depend 
on specific geometry but would not be justified in general, as our result shows. 



Comparison with previous numerical simulations (see, e.g., Giacalone & Jokipii (1999)) 



requires the evaluation of the MFLRW, for isotropic turbulence. Using Eq. (25), we find for 



the average square displacement of MFL (same contribution along x- and y-axes) 



d 



MFL 



{t)k 



1 fSBV q-1 



- 



(46) 



here the term integrated over scales smaller than coherence scale 27r /kmin (or k > kmin) is 
recast as 

,-2 



1 y 

H\{y,q) = — {cosy + ySi{y)) + ——sin y 
qy 2 + q 

g2 + 2g+l 

- 2 — -1 - q, y) 

q{2 + q) 



(47) 
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and the term integrated over scales larger than coherence scale 27r /kmin (or ^ < ^mm) is 
recast as 



yo 



^™ , Si(?/) fSi{y) siny cosy\ 



y 



yo 



(4J 



Figure |4] shows the diffusive behaviour of MFL for a magnetic turbulence with isotropic 
wave-number spectrum. As for the slab case, scales larger than coherence length {k < kmin 
dominate over the turbulent contribution (see Fig. 111). The leading terms in the diffusive 



limit of H2{y, q) are found in the integral and in the first term in parenthesis of Eq. (48). By 



integrating by parts and using Gradshteyn & Ryzhik (1973), Eq. (4.421.1), it can be found 



{t)k 



1 ^6By q-l 



Hy)^og{y)Z 



yo 



K 



kn 



(49) 



where C = 0.577215 is the Euler constant and ym = kminV\\t — kmin^g^t {yo ~ korgflt). 
We find that the MFL of a 3D isotropic weakly turbulent magnetic field are superdiffusive 



according to Eq. (49). This result does not imply the superdiffusion of the particles which 



propagate diffusively in a 3D isotropic turbulence (Giacalone & Jokipii 1999). Before being 



transported superdiffusively along a field line, a particle will undergo parallel scattering, 
not taken into account in this paper, and eventually decorrelate to another field line. These 
results might be qualitatively extended to MHD-turbulence with a small k-anisotropy so 
that results in diffusion regime found here still apply. In this case, particle drift retains its 
diffusive character shown in Fig. [3} even if only scales smaller than the coherence length 
{k > kmini in dotted-dashed) are taken into account. Therefore, at small length-scales, the 
perturbative approach based on 5B <^ Bq could be applied to solar wind turbulence. At 
length-scales much larger than the coherence scale {k <^ kmin), where the approximation 
6B ^ Bq breaks down in the solar wind, these conclusions cannot be extended. 



5. Discussion and conclusion 



We have described analytically the time evolution of individual charged particles drift 
and MFLRW across a static magnetic field to first-order in the magnetic fluctuations. We 
consider the case where the motion perpendicular to the average magnetic field is dominated 
by guiding center motion which includes the meandering of the MFL and the drift from the 
first-order orbit theory, in the approximation that the particle gyroradius is much smaller 
than the length scale of magnetic field variations. In contrast to previous models for the 
perpendicular transport, we do not assume diffusive scattering at all times; this allows us 
to treat consistently the slab turbulence perpendicular diffusion. Drift and MFL transverse 
transport are explicitly computed for both the slab and 3D isotropic cases. In the slab 
case, the time-evolution of the drift displacement shows how the particle diffusion from the 
MFL is suppressed. The instantaneous slab drift coefficient diffusion transverse to the local 
field depends on the turbulence power-law spectral index; for a Kolmogorov turbulence is 
found to decrease as t~^/^, slower than compound diffusion displacement which is 

however computed transversally to the average field and not to the local magnetic field 
as in this paper. The recovery of the MFL coefficient diffusion of QLT shows that this 
result does not depend on assuming MFL parallel diffusion. Secondly, we provide analj^ical 
time- dependence of drift and MFL coefficients of diffusion for a 3D isotropic turbulence. 
We found that for a 3D isotropic turbulence the particle drift from the local field hne is 
diffusive, whereas the field line itself is superdiffusive. Previous numerical simulations are 
not contradicted by our result which is obtained neglecting scattering parallel to the mean 
field. For the slab, we find that MFLRW is dominated by length-scales larger than the 
coherence length whereas particle drift from the local field line is dominated by length-scales 
smaller than the coherence length. This disentaglement does not hold for 3D isotropic 
turbulence: MFLRW is still dominated by large length-scales whereas as far as drift is 
concerned turbulent energy contribute at all scales, both below and above the coherence 
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length. The study carried out here provides a framework for particle transport in the solar 
wind and supernova remnant blast wave turbulence and questions the common assumption 
that cosmic-rays trajectory follow the magnetic field line. 
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Fig. 1. — Average transverse drift coefficient in units of (rescaled by a factor 10^) as a 
function of fit for k'^'^'^rg = lO'^, fc^^'^/fc™" = 10^ A;™"/A;f = 10^, q = 5/3 and 6B/B = 0.1. 
The diffusive approximation in Eq. ( [32| is superposed (dashed line). The departure from the 
perturbed field line of the guiding center, which in the ffist-order orbit theory represents the 
real particle motion with a good approximation, decreases to zero. Thus in a slab turbulence 
charged particles are tied to the weakly perturbed magnetic field hues. 
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Fig. 2. — Average magnetic field line transverse displacement as computed in Eq.(25) is 
compared with the QLT coefficient diffusion Dmfl as a function of Vlt for k^°'^rg = 10^'^, 
l^rnax^j^min ^ ^q4^ A;™«/A;J = 10^, q = 5/3 and 6B/B = 0.1. The horizontal line represents 
the quasi-linear limit. The field-line wandering, by computing the magnetic fluctuation along 
the unperturbed trajectory from times smaller than the time needed to particles to travel 
the turbulence correlation length, reaches asymptotically the QLT limit. 
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Total instantaneous diffusion coefficient 

. k < k^^„ 
Diffusive limit 




10' 10' 

n t 



Fig. 3. — Total average transverse drift is compared with the distinct contributions at length- 
scales smaller and larger than the coherence length in the diffusive regime in units of ft;^ as a 
function oiUt for /c^^^^r^ = IQ-^, kmax/kmin = 10^ kmin/h = 10^, q = 11/3 and SB/B = 0.1. 
For comparison the diffusive limit, k]j/kb, is shown. In contrast to the slab, length-scales 
larger than the coherence length give a non-negligible contribution to diffusion. 
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10 



10° 



Fig. 4. — Total average magnetic field line transverse displacement for a 3D isotropic tur- 
bulence is compared with the distinct contributions at length-scales smaller and larger 
than the coherence length in the diffusive regime as a function of Qt for kmn^Tn = 10^^ 



10'', kmin/ko = 10^, q = 11/3 and 6B/B = 0.1. The diffusive limit, Hmfl^t. 



defined in Eq. (49), is overlaid (cf. text) 



